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This contribution provides a pedagogical introduction to and review of the current status and ongoing progress in the 
development of Monte Carlo tools for the calculation and simulation of high-Q^ processes in hadronic collisions. 



1 Introduction 

In April 2007 the Large Hadron Collider (LHC) will 
start providing pp collisions at = 14 TeV. The 
main goal [ ^ of this new enterprise is the explo- 
ration of a yet uncharted energy domain, and QCD 
will provide an essential tool for the analysis and the 
decoding of the immense set of data that will be col- 
lected by the three detectors currently being built [E]- 
The production rates for most SM particles and pro- 
cesses are mindboggling. Huge statistics of final states 
nowadays totally unaccessible at the current accelera- 
tors will be available, and will allow measurements of 
unprecedented accuracy and depth [El-lEl- Table Q] 
gives few examples of cross-sections for some of the 
most relevant processes. 

Table 1. Benchmark cross-sections for few SM processes 
at the LHC, L = 10"^''cm~^s~^ . In the case of jets and 
photons, we assume a |77j < 1 rapidity cut. 



Process 


cr (nb) = cvts/s 


Jets, Et > 0.1 (2) TeV 

eue 

Z e+e- 
Photons (Et > 60GeV) 
cc, bb 
tt 


10^ (10-*) 
20 
2 
20 

8 X 10^, 5 X 10'°^ 
0.8 



In addition to direct manifestations of new phenom- 
ena, these measurements could ultimately lead to in- 
direct evidence for physics beyond the SM itself. The 
huge lever arm in energy available through the mea- 
surement of high transverse-energy {Et) jets [|4| will 
allow to probe the smallest distance scales ever ac- 
cessed. One year of high-luminosity running will give 
tens of events with jets with Et > 3 TcV. Compared 
to the Tevatron, where jets up to 600 GeV will be 
observed, this is a factor of 5 increase in the scale at 
which the quark form factor can be explored. 

The immense samples of EW gauge bosons will en- 
able high-precision measurements of the W mass 
(±15 MeV) and of the gauge bosons selfcouplings [ 



IS). The study of Drell-Yan final states will be sensi- 
tive to several possible new phenomena. For example, 
possible contact interactions mixing light quarks with 
leptons will be probed up to scales of the order of 25- 
30 TeV, well in the region where new strongly interact- 
ing phenomena related to the EWSB may take place. 
New U{1) gauge bosons with Z-like couplings will be 
observed up to masses of the order of 4-5 TeV. The 
large statistics will allow measurements of the W and 
Z total cross sections with accuracies significantly bet- 
ter than 1%, providing luminosity monitors which are 
only limited by the precision of the theoretical predic- 
tions for these rates (predictions currently estimated 
to be around 5%, dominated by uncertainties in the 
partonic densities of the proton [W)- 

Top quarks will be produced in great abundance, at a 
rate of approximately 1 pair/s. Precise determinations 
of the top decay properties and the top mass with 
an uncertainty of about 1 GeV will be possible [ |B] . 
The study of flavour will be enriched by a thorough b- 
physics programme [|7]. The tens of billions of bottom 
quarks produced [|Hj, will allow to pin down with great 
accuracy and redundancy the CKM matrix elements, 
to probe in full detail the parameters of C^ in the Bq 
systems {q — u,d,s), and to study rare decays with 
branching ratios at the level of 10~^. 

In all these cases, it is required that the best possi- 
ble control on the QCD production and decay mech- 
anisms be available. For this reason, in the recent 
years a great effort has been put into the development 
of tools enabling the description of the final states re- 
sulting from high-energy pp collisions. These tools go 
under the common name of Monte Carlo (MC) codes, 
since the state-of-the-art knowledge about QCD is im- 
plemented using numerical MC techniques. MC de- 
velopment is a very technical topic. This review is 
intended to provide a pedagogic and qualitative intro- 
duction to the main issues, and to some of the most 
relevant ideas and topics which are driving the fast 
development in the field. More accurate discussions, 
and all the necessary details, will have to be looked 
for in the references. 
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2 Overview of the available approaches 

The starting point of all QCD analyses of high-Q^ 
processes is summarized by the factorization theorem, 
most clearly expressed by the following relation: 

X j^k 

-^^}^pllF{X^X;Qf), (1) 

where: X represents a given hadronic final state (FS); 
X is an arbitrary partonic FS; fj{x, Qi) is the number 
density of partons of type j carrying the momentum 
fraction x of the nucleon at a resolution scale (factor- 
ization scale) Qi] <jfk{Qi,Qf) is the partonic cross- 
section for the transition between the initial partonic 
state jk and the final partonic state X, considered 
at resolution scales Qi and Qf for initial and final 
state factorization, respectively^; F{X ^ X; Qf) rep- 
resents a transition function from the partonic final 
state to the given observable X. This may include 
fragmentation functions, hadronization effects, as well 
as the result of experimental cuts or jet definitions. 
The sum over final states X can be thought of as a 
sum over all possible histories leading to the same ob- 
servable configuration. In a similar fashion, the par- 
ton density (PDF) fj{x,Q) can be thought of as the 
sum over all initial-state (IS) evolution histories lead- 
ing to the initial-state parton j with momentum frac- 
tion X. According to the factorization theorem, the 
possible IS and FS histories, and their relative proba- 
bilities, are independent of the hard process, and only 
depend on the flavours of the partons involved and 
on the resolution scales. Once an algorithm is devel- 
oped to describe IS and FS evolution, it can therefore 
be applied to the partons arising from the calculation 
of an arbitrary hard process. Depending on the ex- 
tent to which possible IS and FS histories affect the 
value of X, three different realizations of the factor- 
ization theorem are used: cross-section "evaluators" , 
parton-level event generators, and shower MC event 
generators. We shall now review these in more detail. 

2.1 Cross-section evaluators 

In this case only some component of the FS is sin- 
gled out for the measurement, all the rest being ig- 
nored {i.e. integrated over). One example is the 
Drell-Yan (DY) process, where we typically look at 
the inclusive spectrum (mass, pr rapidity, etc.) of a 
leptonic FS. In this approach, no event needs to be 

^An additional dependence on the renormalization scale /^r is 
also present, but this is usually associated and identified with 
Qi and/or Qf. 



generated, it is enough to define the variables relative 
to the considered object (the lepton(s), or a jet). Ex- 
perimental selection criteria {e.g. a jet definition or a 
detector acceptance) are applied to parton-level quan- 
tities. Provided these are infrared and coUinear finite, 
it does not matter what F is, as we assume that all 
histories of the final partonic state lead to the same 
observable, and integrate to 1 because of unitarity: 
J2x P{^^^) = 1- Thanks to the inclusiveness of the 
result, it is conceptually "straightforward" to include 
higher-order corrections, as well as to resum classes 
of dominant and subdominant logarithms. Next-to- 
leading order (NLO) results are known [ 0] ^] for 
most processes both within and beyond the Standard 
Model (SM). In addition, next-to- next-to-leading or- 
der (NNLO) cross-sections have been calculated for 
the DY-type processes [QT^ and for Higgs production [ 
I12| . For a review of progress in NNLO calculations for 
higher- multiplicity final states, see [ I13j . 

2.2 Parton-level event generators 

Here parton- level (PL) configurations {i.e. states with 
quarks and gluons) are generated, with probabilities 
proportional to the relative perturbative matrix el- 
ement (ME). The transition function between a FS 
parton and the observed object (jet, missing energy, 
lepton, etc.) is unity, so there is no need to model the 
exclusive realization of the histories associated to the 
PDF and to F, as they all lead to the same observ- 
able. Experimentally, this is equivalent to assuming 
a smart jet algorithm (that would associate a jet to 
each hard parton) and linear detector response (the 
energy and direction of a measured jet will not depend 
on its inner structure). The advantage over the cross- 
section evaluators is that, with the explicit representa- 
tion of the kinematics of all hard objects in the event, 
more refined detector analyses can be performed, im- 
plementing complicated cuts and correlations which 
are otherwise hard to simulate with the inclusive ap- 
proach. PL event generators are typically used to de- 
scribe final states with several hard jets. Due to the 
complexity of the ME evaluation for these many-body 
configurations [ I14L 1151 ITC) , calculations are normally 
available only for leading-order (LO) cross-sections. 
In this case, several tools [ 1^ have recently 
become available, covering all of the necessary pro- 
cesses for signal and background LHC studies, with 
jet multiplicities all the way up to 4, 5 or 6, depend- 
ing on the specific process: {W .ff')QQ+ jets 
{Q being a heavy quark, and / — £,q)] {Z/j* — s- 
//) QQ+ jets (/ ^£,iy)]{W~* ff) + charm -f jets 
(/ = £, g); {W ^ //') + jets and (^7*_; //)+ jets; 
nW + mZ + lH+ jets; QQ+ jets; QQQ'Q'+ jets, with 
Q and Q' heavy quarks; HQQ jets; N jets; A^7-|- jets. 
Most of these processes have been evaluated by more 
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than one group, providing additional robustness to 
the calculations and their numerical implementation. 
NLO PL event generators are also available for sev- 
eral low-jet-multiplicity final states, as discussed in 
Sect.EU 



2.3 Shower MC event generators 

Shower MC generators [ provide the most com- 
plete description of the FS. Their goal is to gener- 
ate events consisting of physical, measurable hadrons, 
with a correct description of their multiplicity, kine- 
matics and flavour composition. These final states can 
therefore be processed through a complete detector 
simulation, providing the closest possible emulation of 
a real event. Shower MC codes such as HERWIG [ I26j . 
PYTHIA [EZj, IS A JET [EHI or ARIADNE [EH] have 
been known and used for several years now, and new 
tools, SHERPA, are becoming available [I3U|. 

After the generation of a given PL configuration (typ- 
ically using a LO ME for 2 —> 1 or 2 ^ 2 processes), 
all possible IS and FS histories ("showers") are gen- 
erated, with probabilities defined by the shower algo- 
rithm. By algorithm, we mean a numerical, Markov- 
like evolution, which implements within a given ap- 
proximation scheme the QCD dynamics. This includes 
the probabilities for parton radiation (gluon emission, 
or 5 — > splitting), an infrared cutoff scheme, and a 
hadronization model. The radiation probabilities are 
defined as exclusive quantities, unitarized by the in- 
clusion of Sudakov form factors. This means that the 
shower evolution itself does not alter the overall cross- 
section, as estimated from the ME evaluation for the 
initial hard process. Therefore a shower MC based 
on LO matrix elements cannot provide an estimate of 
the K factor. Radiation probabilities implement lead- 
ing soft and coUinear logarithms, plus some sublead- 
ing classes of logs. Quantum-mechanical correlations 
between different emissions are negligible in the case 
of coUinear emission, since in this case gauges can be 
defined where the interference between different dia- 
grams is numerically suppressed; in the case of soft 
emission at large angle, the association to a specific 
emitter is ambiguous, as the interference between dif- 
ferent emission diagrams is large. Fortunately, soft 
emission at large angle is heavily suppressed, as the 
interference effects between different diagrams lead to 
destructive interference. Quantum coherence can then 
be implemented in a branching evolution via, for ex- 
ample, an angular ordering prescription [ I31|: if Oij 
is the angle between two colour-connected partons, 
the soft-gluon emission probability from the pair i, j is 
given by the sum of two independent emission terms, 
constrained within angular regions close to the two 



emitters: 

p(0,)e(% - Oi) + p(0,)e(% - Oj) (2) 

where 0i, 9j are the angles between the soft gluon and 
the colour connected partons emitting it, and P{9) is 
a positive definite probability. 

3 Complementarity of the three ap- 
proaches 

Each of the above three approaches has its virtues and 
shortcomings. Until recently, we had to select the ap- 
proach that was closest to our needs, and live with the 
fact that it had drawbacks. We first outline a compar- 
ison between the merits of the three approaches, and 
then review the recent progress in merging them. 

3.1 Final-state structure 

Shower MC's provide the most complete description. 
The full information about the event is available, in 
terms of physical hadrons. So this is the only tool that 
allows realistic detector simulations, crucial when it 
comes to measuring with precision some quantity, such 
as a particle mass or a coupling strength. In the case 
of PL MC's only the information on the hard partons 
is given; this is good enough for studies with naive de- 
tector simulations based on pure geometry and exact 
energy resolution. This is an excellent tool to quickly 
study systematics such as the PDF dependence of an 
acceptance, but cannot be used, for example, as a tool 
to reconstruct the top quark mass from real W+4: jet 
events. Cross-section evaluators typically don't pro- 
vide much information on the structure of the final 
state. 

3.2 Higher-order corrections: real, hard emis- 
sions 

In the shower MC's, emissions are treated within the 
coUinear or soft approximation. In particular, the 
angular-ordering prescription, introduced to imple- 
ment quantum coherence effects for large-angle gluon 
radiation, will suppress hard large-angle emissions as 
well. Since large-angle hard emissions are not con- 
strained by angular ordering, the shower approxima- 
tion will typically underestimate the rate of multijet 
final states. On the contrary, PL MC's have all dia- 
gram interference terms among multiple hard partons 
taken into account, and provide the best tool to de- 
scribe multijet final states. 

3.3 Higher-order corrections: virtual effects 

Virtual diagrams are included in the soft approxima- 
tion in the shower MC's. They are needed to cancel 
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the infrared (IR) singularities due to the real emis- 
sion of soft gluons, and appear through the Sudakov 
form factors to enforce the unitarity of the shower 
evolution. Virtual effects are exactly included in all 
NLO cross-section evaluators. In addition to the ex- 
act treatment of real emission, NLO (or NNLO) cross- 
section evaluators provide the most accurate determi- 
nation of inclusive rates. In order to implement them 
in ME event generators, however, one must deal with 
the problem of negative- weight events. These arise 
to enforce the cancellation of the positive-weight in- 
finities present in the soft and coUinear real emission 
diagrams, with negative- weight infinities associated to 
the interference of the virtual and tree-level diagrams. 



4.1 Inclusion of NLO corrections in ME gen- 
erators 

At NLO one needs to account for both real-emission 
diagrams and virtual corrections. Virtual corrections 
to an N-body FS give still an N-body FS, while real 
emissions lead to an (N-l-l)-body FS. An NLO event 
generator must therefore generate both N-body and 
(N-l-l)-body FS's. The two however cannot be treated 
separately, since the contribution of the virtual part to 
the cross-section is equal to minus infinity, while the 
real part has divergencies for coUinear and soft config- 
urations. The structure of the cross-section in these 
regions can be parameterized in an idealized form as 
follows: 



3.4 Resummations 

Resummation [ |22l E El of leading and subleading 
logarithms appearing at all orders of perturbation the- 
ory when largely different energy scales appear in the 
event kinematics (for example small-pr production of 
heavy objects, with M ^ pT, or large-p^ production 
of light objects, pr 3> m-), is possible in the context 
of cross-section evaluators. The resummation of log- 
arithms corresponds to the integration over multiple- 
emission and multi-loop diagrams, and is therefore, by 
its nature, an inclusive calculation. As a result, so far 
no ME generator includes them, although new ana- 
lytical [ 1211 IMl and semi-numerical [ Oil approaches 
have been proposed to provide resummation correc- 
tions to arbitrary multi-parton final states. Shower 
MC's incorporate resummations via the multiple emis- 
sions taking place in the shower, and the correspond- 
ing Sudakov form factors, which account for the multi- 
loops. However, due to the unitarity of the evolution, 
these resummations cannot affect the overall rate of 
a process. The resummation of threshold-like loga- 
rithms (such as those entering in the determination of 
the ti cross-section [EH) affect the rates, and so far are 
only included in cross-section evaluator codes. Resum- 
mation of small- a; effects [03; finally, is implemented 
in some shower MC's[OHlEni! mostly developed for 
HERA studies, and is not yet a standard component 
of LHC tools. 



4 A Monte-Carlo of Everything 

The efforts of the past few years have improved con- 
siderably the flexibility of the above tools, and have 
allowed the construction of MC codes which merge the 
merits of different approaches. The main lines of de- 
velopment are summarised, in a very simplified way, 
in this section. 



dx 



fix) 



(3) 



where by definition of ()_ 



da 
dx— 

dx 



(1-^) 



(4) 



The limit a; — > 1 corresponds to kinematical configura- 
tions where the (N-l-l)-body FS degenerates (via soft 
or coUinear emission) to an N-body FS. In the case of 
coUinear emission, for example, x — cos 9. For heavy 
quark (Q) pair production, x — "^q^q/'^ ^ 1 corre- 
sponds to the soft-emission limit. f{x) is continuous 
in X — 1, and therefore any integral over an arbitary 
range of x, including possibily a; = 1, is finite. This is 
the essence of the virtual/real cancellation of infinities. 
If we formally remove the integration from eq. 0] we 
can "define" a non-singular differential cross-section 
as: 



da 
dx 



/W-/(l) 
(1-x) 



(5) 



One can interpret this relation as follows: for any given 
(N-l-l)-body kinematical configuration C{x) {x ^ 1), 
whose weight is given by /(x)/l — a;, we can associate 
an N-body virtual configuration C(l) with weight 
— /(1)/1 — a;, whose kinematics is obtained by the 
a; 1 limit of C(a;). These virtual configurations are 
typically called "counter-events" . One can construct 
an event generator by adding to each (N-l-l)-body final 
state its corresponding N-body counter-event. Both 
events will have a finite weight, although the counter- 
event's is negative. Since their kinematics is differ- 
ent, event and counter-event will typically populate 
different bins of our histograms. For bin sizes suf- 
ficiently large (namely for observables sufficiently in- 
clusive) , summing events over the full phase-space will 
nevertheless lead to positive rates in each bin. This is 
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not guaranteed to happen when the bins are too smalL 
For example, let us consider the integral of the cross- 
section in a bin covering the range 1 — e < a; < 1, with 
e small: 

The first contribution is from the A^ + l-body FS's, the 
second from the counter-events which all accumulate 
at X = 1. Simple algebra leads to the following result: 

/ da^C + f{l)loge, (7) 
Ji-t 

where C is a finite constant when e — > 0. Since /(I) 
is positive, when the bin-size e is small enough the in- 
tegral becomes negative. This is an indication that 
radiative corrections in this small corner of phase- 
space are large. In other words, if we try to push 
the calculation in a region of low-inclusivity, probing 
the final-state structure with very fine resolution, the 
fixed-order perturbative approximation breaks down. 
Higher-order corrections become very large, and have 
to be included to restore the positivity of the cross- 
section in that bin. The smaller the bin, the larger 
the number of orders required. 

In concrete applications life is complicated by the iden- 
tification of the functions f{x), the interplay of soft 
and coUiner singularities, and the description of the 
phase-space in terms of suitable variables; however 
the above simplified description captures the main fea- 
tures of the technique. It has been used, in various dif- 
ferent implementations, for the development of NLO 
ME event generators covering several of the interest- 
ing LHC processes. Examples include: 2-jets [1301 and 
3-jets [|4]j, heavy quarks [02, and vector boson [021 
production. For a detailed list of available tools, see [ 
HI) The extension of these techniques to NNLO calcu- 
lations has yet to be formulated. Work is in progress [ 
IE], but the difficulty is immense, and it will still be 
some time before we can implement the already known 
NNLO matrix elements for DY and Higgs production 
into an event generator. 

4.2 NLO corrections in shower MC's 

The necessity to include NLO corrections in ME gen- 
erators is twofold. On one side, only shower MC's 
provide a representation of the final state complete 
enough to allow realistic detector simulations. Inclu- 
sion of the NLO matrix elements for the hard pro- 
cess, which will provide cross-sections with full NLO 
accuracy, is a natural improvement of these essential 
tools. On the other hand, as mentioned in the previous 
subsection, the inclusion of NLO effects in fixed-order 



ME MC's leads to distributions which are not posi- 
tive definite, thus calling for a tool where these large 
(and possibly negative) logarithmic effects which arise 
at any fixed order in some corners of phase-space can 
be properly resummed. This goal can be achieved via 
the inclusion of the NLO ME's in the shower MC. 
A priori one may expect this task to be ill defined, 
as shower MC's already incorporate part of the NLO 
effects: they have real emissions, as well as virtual ef- 
fects included in the Sudakov form factors. The naive 
introduction of NLO ME's would then lead to dou- 
ble counting. This is what kept people skeptical for 
many years about the viability of a NLO shower MC. 
Brilliant work by Frixione and Webber recently 
showed how this merging can be done very effectively 
in what they called a MC®NLO. One starts by identify- 
ing the analytic form of the approximation used by the 
shower MC to describe real emission and the leading- 
order virtual correction contained in the Sudakov form 
factor. One can then subtract these expressions from 
the NLO matrix elements; since the shower approx- 
imation has the correct residue for all singular con- 
tributions, the subtracted NLO matrix elements are 
finite. In the simple formalism used in the previous 
section, one can represent the subtraction from the 
real emission term as follows: 

da ^ f{x) - fMcjx) ^g-j 
dx (1 — x) ' 

where fuci^) is the approximate MC expression for 
the real emission matrix element, with the condition 
that fMc{^) = /(!)• In this way the x 1 singularity 
is not removed by the merging with the virtual cor- 
rection, but by letting the shower algorithm handle it 
and absorb it into the Sudakov form factor. As for the 
virtual part, the singular contribution is all contained 
in the shower approximation, and what is left for the 
NLO correction to describe is just a finite term with 
the N-body, Born-like, kinematics. One is still left 
with positive and negative weight events, since the dif- 
ference between the exact non-singular terms from the 
full NLO calculation and those used in the shower can 
have either sign. However, since the residual positive 
and negative weights are bounded, one can define an 
unweighting procedure whereby positive- weight events 
are unweighted against the maximum positive weight, 
and negative- weight events are unweighted against the 
minimum negative weight. This procedure has been 
successfully implemented in MC@NLO codes [ 0^1 QH] 
describing heavy-quark pair, Higgs, DY and gauge bo- 
son pair production. There is no obstacle to extend- 
ing it to all remaining processes known at NLO. Other 
approaches have also been developed into alternaitve 
codes for DY [47^ and vector boson pair produc- 
tion. 
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The inclusion of NLO corrections in the shower MC 
guarantees that total cross-sections generated by the 
MC reproduce those of the NLO ME calculation, 
thereby properly including the K factors and reducing 
the systematic uncertainties induced by renormaliza- 
tion and factorization scale variations. At the same 
time, however, the presence of the higher-order cor- 
rections generated by the shower will improve the de- 
scription of the NLO distributions, leading to depar- 
tures from the parton-level NLO result. This is shown, 
for example, in fig. ^ which shows the pt spectrum 
of a tt pair resulting from the pure NLO calculation, 
from the LO shower, and from the MCQNLO improve- 
ment. At large pxi a region dominated by the NLO 
effects, MC@NLO faithfully reproduces the hard, large- 
angle emission distribution given by the NLO matrix 
elements. At small pT, a region dominated by mul- 
tiple radiation and higher-order effects, the MCQNLO 
departs significantly from the NLO result, while prop- 
erly incorporating the Sudakov resummation effects 
only available via the IS shower evolution. 



.-^ in-l 




1 2 

logio(pr/GeV) 



Figure 1. Transverse momentum distribution of top 
quark pairs using three different approaches: the LO HER- 
WIG MC, the parton-level NLO MC, and the merging of 
the two into MC@NLO. Figure from [HH]. 



4.3 Merging multijet ME generators and 
shower evolution 

The inclusion of NLO ME's in the shower MC's guar- 
antees the correct description of the emission of one 
extra hard parton (ultimately a jet) from the Born, 
LO process. In the case of NLO corrections to the dijet 
final states, for example, this means that all topologies 
with up to three jets will be accurately described. To 
go beyond this in a NLO framework, however, requires 



the knowledge of NLO ME's which are not available 
today, and won't be for a longtime to come in the 
case of the highest multijet FS's which are of inter- 
est for several LHC studies. As mentioned earlier, 
the description of multijets obtained from the shower 
evolution is inaccurate, since hard radiation at large 
angle is suppressed by the angular ordering prescrip- 
tion. One therefore needs an approach in which multi- 
parton events generated using the exact LO ME gener- 
ator can be consistently evolved into multi-jet FS's via 
a shower MC. As in the case of the inclusion of NLO 
corrections, the main problem to be dealt with is that 
of double counting. I shall describe this problem by 
discussing a specific example [ 49 : consider inclusive 
production of 3-jet events, with jets defined by a cone 
of size Rj in 77 — space and transverse energy larger 
than E™"\ We can generate these events by generat- 
ing parton-level configurations with partons separated 
by AR > Rj, and let the shower evolve them into jets. 
By and large, there will be a one-to-one correspon- 
dence between the generated hard partons and the 
jets, and the angular distributions of the three jets will 
include correctly all interference effects among the var- 
ious diagrams. There will be configurations, however, 
where the correspondence is not guaranteed. Take for 
example events where two jets have Et larger than 
the third, and the third itself is well above the i?™™ 
threshold: En Et2 > Ets > EJp'". These events 
can be generated in two independent ways. On one 
side we can start from three partons with the kine- 
matics of the three jets. After evolution, the partons 
will generate the desired jets. On the other, we can 
start from configurations where the two leading par- 
tons have energies of the order of Eti and Et2, but 
the third can be as soft as i?™™. These events will 
typically not have a third jet with Et Ets, since 
the generated parton will be softer. However, hard ra- 
diation by one of the two leading jets may lead, with 
probability as , to the generation of an extra jet with 
the desired energy. In other words, the shower evolu- 
tion in these cases could produce jets with transverse 
energy larger than that of partons already present in 
the hard event. While the probability of this happen- 
ing is parametrically of order as relative to the LO 
process, and could therefore be considered a higher- 
order correction compatible with the LO approxima- 
tion of the ME calculation, configurations with two 
hard partons and a much softer one arc enhanced by 
large logarithms. This is the result of the large phase- 
space available for the emission of a soft parton in a 
hard event. One can estimate that this enhancement is 
of order log(£'T3/-E'™™). As a result the overall prob- 
ability that the third jet be emitted from the shower 
becomes a number of order ag log(£'T3/-£'™'"). On 
one side this could be numerically of order 1, therefore 
turning an NLO effect into a LO one. On the other. 
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this dependence on iJ™™ of the cross-section for 3- 
jet events well above the E"™™ threshold is a totally 
unphysical result! This paradox is caused by the dou- 
ble counting of equivalent configurations which this 
approach gives rise to. To solve the paradox, and re- 
move the E^™ dependence, one has to carefully check 
that a given phase-space configuration is generated 
only once: either directly by the parton-level event, 
or by the shower evolution. 

Some approaches to this problem, aiming at different 
levels of accuracy, have been introduced recently . The 
first [ inni is generically known as "matrix-element 
correction" technique (MEC), as it corrects the ap- 
proximate ME for the emission of the hardest gluon 
in a given process by using the exact LO ME. The 
second is known as CKKW [E]; its goal is to imple- 
ment multi-jet ME corrections at the leading (LL), or 
next-to-leading (NLL) logarithmic level. In the MEC, 
one starts by identifying analitically the phase-space 
region covered by the shower algorithm. In the 
case of e+e^ qqg, for example, this is given by 
a subset of the full phase-space domain A defined by 
A = [1 < xi -ha;2 < 2] n [a;i < 1], where Xi = 2E,/\/S. 
Q. contains the singular regions corresponding to soft 
{xi = 0) and coUinear {xi = 1) gluon emission. The 
integral ctqo of the 3-body cross-section over the com- 
plement of fi, = A — 51, is therefore free of sin- 
gularities and finite. The integral over Q, an, is also 
finite, once the virtual corrections at the edge of phase- 
space are included. The MEC MC generation then 
works by deciding on an event-by-event basis whether 
to generate the event in fl'^ or in fl, based on the 
relative value of the respective cross-sections. If the 
event falls in f2 one generates a LO e+e^ — s- qq event. 
If it falls in 17°, one generates a e+e^ qqg event. 
In both cases, the events are then evolved through 
the shower. Small adjustments should then be made 
in the first case to ensure a proper continuity across 
the boundary between the two domains. This tech- 
nique has been applied also to DY production [ , 
and to top decays [[SSj. Its extension to more compli- 
cated processes, however, is made particularly difficult 
by the need to provide an analytic description of ft. 
When there are more than 3 coloured partons in the 
process (as e.g. in dijet or heavy quark pair produc- 
tion at the LHC), this becomes very hard and imprac- 
tical. The CKKW approach circumvents this problem 
by limiting its precision goal to a LL accuracy (NLL 
for e~^e~ —^ multijets). In this approach the double 
counting is removed not by exactly separating a priori 
the domains fl and fl", but by a probabilistic rejection 
procedure applied to events falling in the overlap, so 
as to ensure that a given phase-space configuration is 
only counted once. One starts by generating samples 
of multi-parton events of different multiplicities, using 



the exact LO ME. The generation is carried out in a 
phase-space domain defined by a Durham-like jet algo- 
rithm suitable for hadronic collisions ; a resolution 
variable for two partons is defined by: 

fcy = inm{ET,i,ETj)Rij , (9) 

if both i and j are in the FS, or by kij = ET,i if J is in 
the IS. Rij is a measure of separation in the transvserse 
plane, for example the standard ^/A^^ + Att^ mea- 
sure. An A^-parton FS is then classified as an iV-jet 
event if fcab > fccut for all possible parton pairings. 
kcut is a resolution threshold, introduced as a param- 
eter necessary to separate the generation of the events 
in samples of different jet multiplicity; the final cross- 
setion, however, should be indepepdent of its specific 
values. Events are extracted from the different iV-jet 
samples with probability proportional to the sample 
cross-section. The jet algorithm can then be used to 
define a tree structure for the event. The two partons 
I, J with the smallest kij are clustered into a single vir- 
tual parton provided the parton types and flavours 
of i and j can be merged. The procedure is repeated 
after removing «, j from and adding i to the list of par- 
tons. The clustering continues until one gets a 2 ^ 1 
or 2 ^ 2 process. The resulting tree is then inter- 
preted as a shower configuration, but with a weight 
given by the exact ME. To be used as an exclusive FS, 
and to allow the successive evolution via the shower, 
the event weight needs to be corrected with the inclu- 
sion of Sudakov form factors for each vertex, and by 
a rescaling of the values of as reflecting the choice of 
scale for made by the shower. This reweigthing 
factor, which can be constructed so as to always be 
smaller than 1, is then used as a probability to keep 
or reject the event. The events kept are then show- 
ered, and shower splittings at a scale larger than kcut 
are vetoed, to avoid the duplication of configurations 
which will otherwise be present in ME samples correp- 
sonding to higher N . In the case of e+e~ collisions [ 
1511 ISni one can prove that this algorithm correctly re- 
produces the weight of an event to NLL accuracy, and 
to rates independent of kcut- In hadron collisions [ 
I56| such a proof is still missing, but this framework 
provides a very good starting point for further devel- 
opments. The results of the first studies in the case of 
M^-l-jets production can be found in ['Ff. 

An additional improvement in the description of mul- 
tijet final states may soon come from the development 
of new shower algorithms, being developed in the con- 
text of the new generation of C-I--I- codes, such as HER- 
WIG-h-h [EH and SHERPA [EOl- Preliminary results 
for e+e~ collisions, comparisons with LEP/SLC data 
and examples of the improved description of multijet 
final states, can be found in [ :58.. Similar work is 
under way for the C-|— I- version of PYTHIA ["SIT. 
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5 Conclusions 

I hope to have suceeded in giving a flavour of the big 
progress and flourishing of new ideas which is cur- 
rently taking place in the field of MC development for 
hadronic collisions. Tools only dreamt of few years ago 
are now available, encouraging and justifying greater 
ambitions. Progress is also taking place in other di- 
rections which I have not had time to cover, such as 
the progress in the evaluation of NNLO matrix ele- 
ments [El J in the description of power corrections [ 
l()l| , underlying event and multiple interactions [ lf)2j , 
hadronization [ I63| , as well as on the uniformization 
of input/output formats for the merging of PL and 
shower generators [ El] The validation of MC tools 
against Tevatron data is also an area of active and suc- 
cessfull research [ EHl EE] ■ Finally, great progress has 
been achieved in the determination of parton densities 
and their sistematic uncertanties , an essential in- 
put for all event simulations. Hoping that the progress 
will continue, and that all of these tools will receive 
proper validation using the results available from the 
Tevatron, I look forward to their testing at the LHC. 
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